QCD with two flavors of Wilson fermions: 
The QCD vacuum, the Aoki vacuum and other vacua 
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We discuss the vacuum structure of QCD with two flavors of Wilson fermions. We derive two possible sce- 
narios: (i) If the spectral density pj/(A., k) of the overlap hamiltonian in a fixed background gauge field is not 
symmetric in X, Hermiticity is violated and Hermiticity violation effects could influence numerical determi- 
nations of the r| meson mass if we are not near enough to the continuum limit, where Hermiticity should be 
recovered; (ii) otherwise we argue that, under certain assumptions, new phases appear beside the Aoki phase, 
which can be characterized by a nonvanishing vacuum expectation value of i\(/,,Y5V|/u + i^idls^d^ ''nd with vac- 
uum states that cannot be connected with the Aoki vacua by parity-flavor symmetry transformations. Quenched 
numerical simulations suggest that the second scenario is more likely realized. 

PACS numbers: 11.15.Ha 11.30.Hv 11.30.Qc 



I. INTRODUCTION 

Since the first numerical investigations of four-dimensional 
non-abelian gauge theories with dynamical Wilson fermions 
were performed in the early 80's 1 1, 2], the understanding of 
the phase and vacuum structure of lattice QCD with Wilson 
fermions at non-zero lattice spacing, and of the way in which 
chiral symmetry is recovered in the continuum limit, has been 
a goal of lattice field theorists. The complexity of the phase 
structure of this model was known a long time ago. The ex- 
istence of a phase with parity and flavor symmetry breaking 
was conjectured for this model by Aoki in the middle 80's 
ijSi, ^ . From that time on, much work has been done in order 
to confirm this conjecture, to establish a quantitative phase di- 
agram for lattice QCD with Wilson fermions, and to delimit 
the parameters region where numerical calculations of phys- 
ical quantities should be performed. References i3ll- i20ll are 
a representative but incomplete list of the work done on this 
subject. 

In this paper we analyze the vacuum structure of lattice 
QCD with Wilson fermions at non-zero lattice spacing, with 
the help of the probability distribution function of the parity- 
flavor fermion bilinear order parameters ll2lll ll22ll . We find 
that if the spectral density Pj/(A,,k) of the overlap hamilto- 
nian, or Hermitian Dirac- Wilson operator, in a fixed back- 
ground gauge field U is not symmetric in X in the thermo- 
dynamical limit for the relevant gauge configurations, Her- 
miticity of /v(/„y5V(/„ + /\(/f/y5\(/f/ win be violated at finite p. As- 
suming that the Aoki phase ends at finite p, this would imply 
the lost of any physical interpretation of this phase in terms 
of particle excitations. In addition, the lost of Hermiticity of 
'V»Y5¥" + iWdJsWd finite p also suggests that a reliable de- 
termination of the r| mass would require to be near enough 
the continuum limit. If on the contrary the spectral density 
Pj/(A,,k) of the Hermitian Dirac-Wilson operator in a fixed 
background gauge field U, is symmetric in X, we explain how, 
under certain assumptions, the existence of the Aoki phase im- 
plies also the appearance of other phases, in the same param- 



eters region, which can be characterized by a non-vanishing 
vacuum expectation value of /v(/„y5\(/„ + i^idJsWd, and with 
vacuum states that can not be connected with the Aoki vacua 
by parity-flavor symmetry transformations. 

The outline of this paper is as follows. In Sec. II we de- 
rive the p.d.f. of /\j/y5T3\|/ and v^y^^ and analyze the con- 
ditions to have an Aoki phase with spontaneous parity-flavor 
symmetry breaking. Section III contains our derivation of the 
p.d.f. of the same fermion bilinears of Sec. II but in pres- 
ence of a twisted mass term in the action. We discuss also in 
this Sec. how if the spectral density p;7(X,K) of the Hermi- 
tian Dirac-Wilson operator in a fixed background gauge field 
U is not symmetric in A,, a violation of Hermiticity manifests 
in a negative vacuum expectation value for the square of the 
Hermitian operator /vj/ys^l/ in the infinite lattice limit. In Sec. 
IV we discuss the other possible scenario i.e., we assume a 
symmetric spectral density p;7(X,K) and then derive, assum- 
ing that the Aoki vacuum does exists, that new vacuum states 
should appear. These new vacua, which are not connected 
with the Aoki vacua by parity-flavor symmetry transforma- 
tions, are analyzed in Sec. V. Section VI contains some nu- 
merical results obtained by diagonalizing quenched configu- 
rations in A"^ ,6"^ lattices. The goal of this analysis was to 
distinguish between the two possible scenarios. Section VII 
contains our conclusions. 



II. THE PROBABILITY DISTRIBUTION FUNCTION IN 
THE GIBBS STATE 

The model we are interested in is QCD with two degenerate 
Wilson quarks. The fermionic part of the Euclidean action is 



Sp = \(/1V(k)i|/ 



(1) 



where W(k) is the Dirac-Wilson operator, K is the hopping 
parameter, which is related to the bare fermion mass mo by 
K = 1/(8 +2wo), and flavor indices are implicit (W(k) is 
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a two-block diagonal matrix). The standard wisdom on the 
phase diagram of this model in the gauge coupling p,K plane 
is the one shown in Fig. 1 . The two different regions observed 
in this phase diagram, A and B, can be characterized as fol- 
lows; in region A parity and flavor symmetries are realized in 
the vacuum, which is supposed to be unique. The Gibbs state 
is then very simple in this region, and continuum QCD should 
be obtained by taking the ^ 0, K ^ 1 /8 limit from within 
region A. We will call region A the QCD region. In region 
B, on the contrary, parity and flavor symmetries are sponta- 
neously broken, there are many degenerate vacua connected 
by parity-flavor transformations in this region, and the Gibbs 
state is therefore made up from many equilibrium states. In 
what follows we will call region B as the Aoki region. 
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FIG. 1: Aoki (B) and physical (A) region in the (P, k) plane. Adapted 
from fl^ with courtesy of the authors. 



The order parameters to distinguish the Aoki region from 
the QCD region are /\j>Y5Tj\(/ and /\j>y5\(/, with ij the three Pauli 
matrices. For the more standard election, /' = 3, they can be 
written in function of the up and down quark fields as follows 



«W5V = «V«Y5¥» + '¥</Y5¥rf 
The Aoki phases are characterized by 

(i¥Y5T:3¥) ^ 
(/¥Y5¥) = 



(2) 



(3) 



The first of the two condensates breaks both parity and fla- 
vor symmetries. The non-vanishing vacuum expectation value 
of this condensate signals the spontaneous breaking of the 
SU{2) flavor symmetry down to U{1), with two Goldstone 
pions. Notwithstanding parity is spontaneously broken if the 
first equation in (|3]l holds, the vacuum expectation value of 



the flavor diagonal condensate (/^YsV) vanishes because it is 
also order parameter for a discrete symmetry, composition of 
parity and discrete flavor rotations, which is assumed to be 
realized 

Following the lines developed in ll2lll . we wish to write 
the p.d.f. of the two fermion bilinear order parameters (|2]i. 
The motivation to develop this formalism was precisely the 
study of the vacuum invariance (non-invariance), in quantum 
theories regularized on a space time lattice, under symmetry 
transformations which, as chiral, flavor or baryon symmetries, 
involve fermion fields. Notwithstanding that Grassmann vari- 
ables cannot be simulated in a computer, it was shown in |2l]] 
that an analysis of spontaneous symmetry breaking based on 
the use of the p.d.f. of fermion local operators can also be 
done in QFT with fermion degrees of freedom. The starting 
point is to choose an order parameter for the desired sym- 
metry (9(\(/,¥) (typically a fermion bilinear) and characterize 
each vacuum state a by the expectation value Ca of the order 
parameter in the a state 



Ca = — / {0{x))ad x 



(4) 



The p.d.f. P{c) of the order parameter will be given by 

P{c)=Y,WaHc-Ca) (5) 



which can also be written as 121 



Pic) = /8(i I Oix)d\ 



-c) 



(6) 



the mean value computed with the integration measure of the 
path integral formulation of the Quantum Theory. 

The Fourier transform P{q) / e''''^'P{c)dc can be written, 
for the model we are interested in, as 



Piq) 



J [dU] [d^td\\t] exp{ -Sym + W (k)\(/ +'^Jd'^x 0{x) } 
J[dU] [d^td\\t] exp{-SYM + ¥^(k)\|/} 



(7) 



In the particular case in which (9 is a fermion bilinear of \\f 
and \|/ 



0{x) = \\f{x)d\\f{x) 



(8) 



with any matrix with Dirac, color and flavor indices, equa- 
tion (|7]i becomes 



^ J [dU] [diiidxv] exp{ -Sym + ¥(W (k) + f q)^} 
^^'^ J[dU][d^d\\f]exp{-SYM + Wi^)'^} 

Integrating out the fermion fields in (|9|l one gets 
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P{q) 



j[dU]e-SyMdtiW{K) 
which can also be expressed as the following mean value 



(10) 



/ det(W(K) + ao) ' 
^^^^ ^ \ detW(K) 



(11) 



computed in the effective gauge theory with the integration 
measure 



Notice that zero modes of the Dirac-Wilson operator, which 
would produce a singularity of the operator in ( fTTT i. are sup- 
pressed by the fermion determinant in the integration measure 
(zero mode configurations have on the other hand vanishing 
measure). 

The particular form expected for the p.d.f. P{c), P{q), de- 
pends on the realization of the corresponding symmetry in the 
vacuum. A symmetric vacuum will give 



P{c)^b{c) 



(12) 



n 



(15) 



where V is the number of degrees of freedom (including color 
and Dirac but not flavor d.o.f.) and "kj are the real eigenvalues 
of the Hermitian Dirac-Wilson operator IV (k) — '^^W{k). No- 
tice that whereas Pj,{q) has not a definite sign, PQ{q) is always 
positive definite. 

The q-derivatives of P{q) give us the moments of the distri- 
bution P{c). In particular we have 



d^Piq) 



dq" 



(16) 



q=0 



Since c is order parameter for the symmetries of the lattice 
action and we integrate over all the Gibbs state, the first mo- 
ment will vanish always, independently of the realization of 
the symmetries. The first non-vanishing moment, if the sym- 
metry is spontaneously broken, will be the second. Thus, for 
the particular case 



P{q) = 1 

whereas, if we have for instance a t/(l) symmetry which is 
spontaneously broken, the expected values for P{c) and P{q) 
are l|2lll 

7^(c) = [7t(c2-c2)l/2]-l -C0<C<C0 

Otherwise P{c) —0 

271 J-7C 

the last being the well known zeroth order Bessel function of 
the first kind, Jo{qco). 

In the simpler case in which a discrete Z(2) symmetry is 
spontaneously broken, the expected form is 



CO - -^/\i/(x)Y5V|/W 

X 

C3 = :^1^/¥(x)Y5T3¥W (17) 

X 

we get 




In the QCD region flavor symmetry is realized. The p.d.f. 
of C3 will be then b{cj,) and (cj) = 0. We get then 



P{c)^^^c-co) + ^d{c + co) 

P{q) = cos{qco) (14) 

or a sum of symmetric delta functions {P{c)) and a sum of 
cosines (P{q)) if there is an extra vacuum degeneracy. 

If we can Pj{q),Po{q) the p.d.f. of i^fjsZj^ii and /\j>Y5\(/ in 
momentum space, we have 




which should vanish since parity is also realized in this region. 
Furthermore a negative value of (cq) would violate Hermitic- 
ity of /vj/ysV- We will come back to this point in the next 
section. 

In the Aoki region |4] there are vacuum states in which the 
condensate C3 SvH takes a non-vanishing vacuum expectation 
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value. This implies that the p.d.f. P^c^) will not be 5(c3) and 
therefore (c^) ( fTSl ) will not vanish. Indeed expression (fTSl ) for 
(C3) seems to be consistent with the Banks and Casher for- 
mula (2^1 which relates the spectral density of the Hermitian 
Dirac-Wilson operator at the origin with the vacuum expecta- 
tion value of C3 ifl^ . 

If, on the other hand, (/\(/Y5X|/} = in one of the Aoki vacua, 
as conjectured in |4], (/\j/yg\(/} = in all the other vacua which 
are connected with the standard Aoki vacuum by a parity- 
flavor transformation, since /Vj/YsV is invariant under flavor 
transformations and change sign under parity. Therefore if 
we assume that these are all the degenerate vacua, we con- 
clude that P{cq) — 5(co) and (cq) ~ 0, which would imply the 
following non-trivial relation 



(19) 




y'^o) - 772 [Ltzj^ 



and 



C3 




(22) 



(23) 



Equation ( |23] l is well known. If we take the infinite volume 
limit and then the m limit we get the Banks and Casher 
result 



III. QCD WITH A TWISTED MASS TERM: THE NON 
SYMMETRIC CASE 

In this section we will consider lattice QCD with Wilson 
fermions with the standard action of previous section, plus a 
source term 



(C3)=27ip(0) 



(24) 



J^im,\\f{x)j5l3\\i{x 



(20) 



that explicitly breaks flavor and parity. The flavor symmetry 
is thus broken from SU{2) to U{1). This is the standard way 
to analyze spontaneous symmetry breaking. First one takes 
the thermodynamic limit and then the vanishing source term 
limit. 

We can calculate again the p.d.f. Po{q) and P^iq) of i'^Js^li 
and i\j/Y5T3\(/ with this modified action. Simple algebra give 
us the following expressions 



Poiq) 




2q ■ 

-finit 



(21) 



where again Xj are the real eigenvalues of the Hermitian 
Dirac-Wilson operator and the mean values are computed now 
with the integration measure of the lattice QCD action modi- 
fied with the symmetry breaking source term. 

By taking the q-derivatives at the origin of ^nd P-i{q) 
we obtain 



(co) 




which relates a non-vanishing spectral mean density of the 
Hermitian Wilson operator at the origin with the spontaneous 
breaking of parity and flavor symmetries. 

The first equation in ( l22l i is actually unpleasant since it pre- 
dicts an imaginary number for the vacuum expectation value 
of a Hermitian operator. However it is easy to see that (cq) — 
because it is order parameter for a symmetry of the modified 
lattice action, the composition of parity with discrete flavor 
rotations around the x or y axis. 

Concerning the second equation in ( l22l i. one can see that 
the first and second contributions to (cq) vanish in the infinite 
volume limit for every non-vanishing value of m. The third 
contribution however, which is negative, will vanish only if 
the spectral density of eigenvalues of the Hermitian Wilson 
operator P£/(A-) for any background gauge field U is an even 
function of X. This is actually not true at finite values of V , 
and some authors fs', suggest that the symmetry of the 
eigenvalues will be recovered not in the thermodynamic limit, 
but only in the zero lattice spacing or continuum limit. If we 
take this last statement as true, we should conclude: 

i. The Aoki phase, which seems not to be connected with 
the critical continuum limit point (g^ = o,K = 1/8) lfT6l[T7ll is 
unphysical since the (cq) would be negative in this phase and 
this result breaks Hermiticity. 

ii. In the standard QCD phase, where parity and flavor sym- 
metries are realized in the vacuum, we should have however 
negative values for the vacuum expectation value of the square 
of the Hermitian operator /lij/ysVI/, except very near to the con- 
tinuum limit. Since this operator is related to the r|-meson, 
one can expect in such a case important finite lattice spacing 
effects in the numerical determinations of the r|-meson mass. 

This is the first of the two possible scenarios mentioned in 
the first section of this article. In the next section we will as- 
sume a symmetric spectral density of eigenvalues of the Her- 
mitian Wilson operator Pf/ (X) for any background gauge field 
U in the thermodynamic limit, and will derive the second see- 
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IV. QCD WITH A TWISTED MASS TERM: SYMMETRIC 
SPECTRAL DENSITY OF EIGENVALUES 

In this section we will assume that the spectral density of 
eigenvalues of the Hermitian Wilson operator p;7(A,) for any 
background gauge field U is an even function of A,. In such a 
case equation (|22] | will give a vanishing value for (cq) at any 
value of nit 

{cl)^0 (25) 
Therefore the p.d. f. of cq is 5(co) and 

(/W5¥) = (26) 

for any value of m,, and also in the m, limit. Thus we can 
confirm that under the assumed condition, (/\(/y5i|/) ~ in the 
Aoki vacuum selected by the external source (l20l i. as stated 
in 101; but since /\j/Y5\|/ is flavor invariant and change sign un- 
der parity, we can conclude that (/\(/Y5i|/) — 0, not only in the 
vacuum selected by the external source ( l20l i. but also in all 
the Aoki vacua which can be obtained from the previous one 
by parity-flavor transformations. In order to see the fact that, 
if there is an Aoki phase with parity-flavor symmetry spon- 
taneously broken, the previous vacua are not all the possible 
vacua, we will assume that is false and will get a contradic- 
tion. 

If all the vacua are the one selected by the twisted mass 
term and those obtained from it by parity-flavor transforma- 
tions, the spectral density of the Hermitian Wilson operator 
will be always an even function of X, since the eigenvalues of 
this operator change sign under parity and are invariant under 
flavor transformations. Then the symmetry of P{/(A,) will be 
realized also at nit = in the Gibbs state. Now let us come 
back to expression ( fTSl l which give us the vacuum expectation 
values of the square of i^ys^ and /iij/y5T3\(/ as a function of the 
spectrum of the Hermitian Wilson operator, but averaged over 
all the Gibbs state (without the external symmetry breaking 
source (|20li). By subtracting the two equations in ( fTSl l we get 

This equation would naively vanish, if the spectral density 
of eigenvalues of the Hermitian Wilson operator were an even 
function of X. Therefore we would reach the following con- 
clusion for the Gibbs state 

(^^3) = {cD (28) 

Nevertheless, S. Sharpe put into evidence in a private com- 
munication (developed deeply in [24]) an aspect that we, 
somewhat, overlooked: A sub-leading contribution to the 
spectral density may affect ( |28] | in the Gibbs state (e-regime. 



in xPT terminology), in such a way that, not only (cq), but ev- 
ery even moment of v^y^^ would vanish, restoring the stan- 
dard Aoki picture. The thesis of Sharpe, although possible, 
would enforce an infinite series of sum rules to be complied, 
similar to those found by Leutwyler and Smilga in the con- 
tinuum [25]. We agree that such a possibiUty is open, at least 
from a purely mathematical point of view: In fact sub-leading 
contributions to the spectral density may exist, and conspire 
to enforce the vanishing of all the even moments of the p.d.f. 
of /vj/ys V- However we believe such possibility not to be very 
realistic, and indeed we have physical arguments, which will 
be the basis for subsequent work on the topic, suggesting that 
the Aoki scenario is incomplete. Therefore, we will reason- 
ably assume in the following (|28] ) to be true, in the case of a 
symmetric pt/(A-), assumption that leads us to the conclusion 
that the Chiral Perturbation Theory may be incomplete, for 
the new vacua derived from (l28T l are not predicted in xVY. 

If as conjectured by Aoki and verified by numerical simula- 
tions, a phase with a non-vanishing vacuum expectation value 
of /ij>Y5T3i|/ does exists, the mean value in the Gibbs state (c^) 
inside this phase will be non-zero. Then equation (l28T l tell us 
that also (cq) will be non-zero inside this phase. But since in 
the Aoki vacua (cq) = 0, this is in contradiction with the as- 
sumption that the Aoki vacua are all possible vacua. This is 
the second possible scenario mentioned in the Introduction of 
this article. 



V. THE NEW VACUA 

To understand the physical properties of these new vacuum 
states we will assume, inspired by the numerical results re- 
ported in the next section, that the spectral density of eigen- 
values Pf/ (A-) is an even function of A, in the Gibbs state of the 
Aoki region (m, = 0). Then equation ( |28] ) holds (taking into 
account the aforementioned discussion raised by S. Sharpe), 
and hence the p.d.f. of the flavor singlet iVj/Ys V order param- 
eter can not be 5(co) inside the Aoki phase, and therefore new 
vacuum states characterized by a non-vanishing vacuum ex- 
pectation value of /\(/Y5\|/ should appear. These new vacua can 
not be connected, by mean of parity-flavor transformations, to 
the Aoki vacua, as previously discussed. 

In order to better characterize these new vacua, we have 
added to the lattice QCD action the source term 

;mf\|>Y5T3i|/ + /eij>Y5V|/ (29) 

which breaks more symmetries than (|20] |. but still preserves 
the U{\) subgroup of the SU{2) flavor. By computing again 
the first moment of the p.d.f. of /vj/ys^l/ and /\(/y5T3\(/ and taking 
into account that the mean value of the first of these operators 
is an odd function of 6 whereas the second one is an even 
function of 9, we get 

,._ , 29/^ -Xj + mj~e^ \ 
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(30) 



where Xj are again the eigenvalues of the Hermitian Wilson 
operator and the mean values are computed using the full inte- 
gration measure of lattice QCD with the extra external source 
( |29] l. This integration measure is not positive definite due to 
the presence of the /vj/ys^l/ term in the action. 

By choosing 9 — nn, in the action and taking the thermody- 
namic limit we get for the two order parameters the following 
expressions 



2rm,A-2 — 2rm^{\ — r^) 
(m2(l-r2) + A.2)^ + 4r2m2A.2 



p{X)dX 



(j\|ry5T3 



2m^{l-r^)+2mtX'^ 
(m2(l-r2)+A.2)^+4r2m2A.2 



p{X)dX (31) 



where p (X) is the mean spectral density of the Hermitian Wil- 
son operator averaged with the full integration measure. 

Taking now the m, limit i.e., approaching the vanishing 
external source ( |29] l point in the 9,OT( plane on a line crossing 
the origin and with slope r, we obtain 



(/iifysv} =2p(0) 



(/\j/Y5T3\|/) -2p(0) 



rt2-r(l -r2) 
(l-r2 + f2)2+4^2j2- 



-dt 



1 -r2 + f2 



(l-r2 + f2)2+ 4^2,2 



dt (32) 



In the particular case of r = (9 = 0) we get the Banks and 
Casher formula 



(/\(/y5T3¥) =27ip(0) 



(33) 



We see how, if p(0) does not vanish, we can get many vacua 
characterized by a non-vanishing value of the two order pa- 
rameters /tj/y5V(/ and /\j/Y5T3X|/. We should point out that the 
value of p (0) could depend on the slope r of the straight line 
along which we approach the origin in the Q,m plane, and 
therefore, even if results of numerical simulations suggest that 
p(0) ^0 when we approach the origin along the line of van- 
ishing slope, this does not guarantee that the same holds for 
other slopes. However the discussion in the first half of this 
section tell us that if p(0) ^ at r = 0, p(0) should be non- 
vanishing for other values of r. 



VI. QUENCHED NUMERICAL SIMULATIONS 

In order to distinguish what of the two possible scenarios 
derived in the previous sections is realized, we have performed 
quenched simulations of lattice QCD with Wilson fermions 



in 4^,6^ and 8^ lattices. We have generated an ensemble of 
well uncorrected configurations for each volume and then a 
complete diagonalization of the Hermitian Wilson matrix, for 
each configuration, gives us the respective eigenvalues. We 
measured the volume dependence of the asymmetries in the 
eigenvalue distribution of the Hermitian Wilson operator, both 
inside and outside the Aoki phase. 

We want to notice that because of kinematic reasons (prop- 
erties of the Dirac matrices), the trace of all odd p-powers of 
the Hermitian Wilson operator W{k) = JsWIk) vanish until 
p — 1, this included. This means that the asymmetries in the 
eigenvalue distribution of W{k) start to manifest with a non- 
vanishing value of the ninth moment of the distribution. We 
have found that these asymmetries, even if small, are clearly 
visible in the numerical simulations. 

In Figs, from 2 to 6 we plot the quenched mean value 



A(P,K,m,) 




(34) 



multiplied by the volume for three different volumes, in order 
to see the scaling of the asymmetries in the eigenvalue distri- 
bution of W(k). As previously discussed,A(P,K,mf) give us a 
quantitative measure of these asymmetries. We have added an 
extra V factor to make the plots for the three different volumes 
distinguishable. Since we found that the value of A(P,k,ot,) 
decreased as the volume increased, the plots of the larger vol- 
umes were negligible with respect to the plot of the smaller 
volume 4^. Multiplying by V all the plots are of the same 
magnitude order 

The nit term in the denominator of ( |34] | acts also as a regula- 
tor in the quenched approximation, where configurations with 
zero or near-zero modes are not suppressed by the fermion 
determinant. This is very likely the origin of the large fluctu- 
ations observed in the numerical measurements of (l34b near 
m, = in the quenched case. That is why our plots are cut 
below m, ~ 0.05; in the physical phase, this cutoff is not re- 
ally needed, but in the Aoki phase it is more likely to find zero 
modes which spoil the distribution. 

Figs. 2 and 3 contain our numerical results in 4^,6^* and S'^ 
lattices at p = 0.001, K = 0.17 and P = 5.0,K = 0.15. These 
first two points are outside the Aoki phase, the first one in the 
strong coupling region. The second one intends to be a point 
where typically QCD simulations are performed. 

Figs. 4, 5 and 6 represent our numerical results in 4^,6''^ 
and 8"^ lattices at P = 0.001, k = 0.30, p = 3.0,k = 0.30 and 
P = 4.0, K = 0.24. These points are well inside the Aoki phase, 
and the structure of the distribution is different from the struc- 
ture observed in the previous plots of the physical phase. Nev- 
ertheless, the qualitative behaviour as the volume increases is 
the same. 

We observe large fluctuations in the plotted quantity near 
m, = 0, specially inside the Aoki phase. However the behav- 
ior with the lattice volume may suggests a vanishing value of 
A(p,K,m() in the infinite volume limit in both regions, inside 
and outside the Aoki phase. If this is actually the case in the 
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FIG. 2: Point outside of the Aoki phase (p = 0.001, K = 0.17) and 
in the strong coupling regime. The superposition of plots clearly 
states that the asymmetry of the eigenvalue distribution decreases as 
^. Statistics: 240 configurations (4"*), 2998 conf. (6"^) and 806 conf. 
(8^) 



FIG. 4: Point inside the Aoki phase ((3 = 0.001, K = 0.30) and in 
the strong coupling regime. Although there is no clear superposition 
of plots, it is evident that the asymmetry goes to zero as the volume 
mcreases. Statistics: 368 conf. (4"*), 1579 conf. (6^*) and 490 conf. 
(84) 




0.02! 



0.02/t 



0.015 - 



0.01 



0.005 




FIG. 3: Another point outside of the Aoki phase (P = 5.0, K = 0.15) 
in a region in which QCD simulations are commonly performed. The 
conclusion is the same as in Fig. 2. Statistics: 400 conf. (4^), 900 
conf. (6"*) and 200 conf. (8^*) 



FIG. 5: Point inside the Aoki phase (P = 3.0, K = 0.30). The asym- 
metry disappears as the volume increases. Statistics: 400 conf. (4^), 
1174 conf. (6"*) and 107 conf. (8"*) 



unquenched model, the second scenario discussed in this arti- 
cle would eventually be realized. 



VII. CONCLUSIONS 

We have analyzed the vacuum structure of lattice QCD 
with two degenerate Wilson flavors at non-zero lattice spac- 
ing, with the help of the probability distribution function of 
the parity-flavor fermion bilinear order parameters. From this 
analysis two possible scenarios emerge. 

In the first scenario we assume a spectral density pu {X, k) 
of the Hermitian Dirac- Wilson operator, in a fixed background 
gauge field U, not symmetric in X. This property is realized at 
finite V for the single gauge configurations, even if a symmet- 
ric distribution of eigenvalues is recovered if we average over 
parity conjugate configurations. We find that under such an as- 
sumption, Hermiticity of /\j/i,y5\(/„ + i^idJiWd will be violated 
at finite p. This lost of Hermiticity for the pseudoscalar flavor 
singlet operator suggests that a reliable determination of the 



r| mass would require to be near enough the continuum limit 
where the symmetry of the spectral density pu (k, k) should be 
recovered. Furthermore assuming that the Aoki phase ends at 
finite p, the violation of Hermiticity obtained in this scenario 
implies the lost of any physical interpretation of this phase in 
terms of particle excitations. 

In the second scenario a symmetric spectral density 
P(/(A,,k) of the Hermitian Dirac-Wilson operator in the in- 
finite volume limit is assumed, and then we show that the 
existence of the Aoki phase implies also the appearance of 
other phases, in the same parameters region, which can be 
characterized by a non-vanishing vacuum expectation value 
of /v(/„y5X|/j, + i'^d'Ys'^ch and with vacuum states that can not 
be connected with the Aoki vacua by parity-flavor symme- 
try transformations. These phases, however, are not related 
to those mentioned in |26], for we keep the twisted mass pa- 
rameter nis equal to zero, whereas the phases studied by G. 
Miinster appear at large m^. 

Sharpe and Singleton ifTolfTsIl performed an analysis of lat- 
tice QCD with two flavors of Wilson fermions near the con- 
tinuum limit, by mean of the chiral effective Lagrangians. In 
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771, 



FIG. 6: Point inside the Aoki phase ((3 = 4.0, K = 0.24). Same con- 
clusions as in the other Aoki plots. Statistics: 398 conf. (4^), 1539 
conf. (6"*) and 247 conf. (8^) 

their analysis, they found essentially two possible realizations, 
depending on the sign of a coefficient C2, which appears in the 
potential energy, expanded up to second order in the quark 
mass term. If ci is positive, a phase with spontaneous flavor 
symmetry breaking and an Aoki vacuum can be identified. If, 
on the contrary, ci < 0, flavor symmetry is realized in the vac- 
uum and a first order transition should appear. Either case may 
be realized in different regions of parameter space. Indeed 
numerical simulations with dynamical fermions, performed 
at small lattice spacing ||27i1, give evidences of metastability 
that can be related with the existence of a first order phase 
transition and hence C2 < 0, while the Aoki phase found at 
smaller P values 1 16] supports C2 > 0. We want to emphasize 
that the new vacua we find are coexisting with the standard 
Aoki vacuum. These new vacua do exist if, and only if, the 
Aoki vacuum exists. But these new vacua are not explained in 
any way in the chiral effective Lagrangian approximation of 
Sharpe and Singleton. On the other hand, in order to recover 
the standard Aoki picture via xPT, an infinite set of sum rules 
for the eigenvalues of the Hermitian Wilson operator must be 
imposed. As we stated previously, this possibility does not 
appeal us, in the sense that it seems unphysical (we would 



call it a 'mathematical' possibility). Nevertheless we have not 
definite proof of the new vacua. 

From our position, these conclusions lead us to cast a doubt 
into the completeness of the Chiral Perturbation theory. In any 
case we believe that, in order to definitely clarify this issue, a 
careful investigation of the spectral properties of Hermitian 
Dirac-Wilson operator for actual gauge field configuration in 
the full unquenched theory is mandatory. 

In order to distinguish what of the two possible scenarios, 
as mentioned at the beginning of this section, is realized, we 
performed quenched simulations of lattice QCD with Wilson 
fermions in 4^,6^* and 8^ lattices, and measured the volume 
dependence of the asymmetries in the eigenvalue distribution 
of the Hermitian Wilson operator, both inside and outside the 
Aoki phase. To measure this asymmetries we diagonahzed ex- 
actly the Hermitian Wilson operator for all the gauge config- 
urations generated with the quenched measure, and measured 
the quenched average A(p,K,m) ( l34b . Due to the fact that con- 
figurations with zero or near-zero modes are not suppressed 
by the fermion determinant in the quenched case, this quan- 
tity fluctuates violently near m = 0, specially for the points 
in the Aoki phase. Notwithstanding that, the observed be- 
havior with the lattice volume seems to suggest a vanishing 
value of A(p,K.m) in the infinite volume limit, both inside 
and outside the Aoki phase. If this were actually the case in 
the unquenched model, the second scenario discussed in this 
article would be realized. However a verification of these re- 
sults in the unquenched case would be very relevant in order 
to discard the first scenario. 
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